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A phase diagram for a quasi-2D metal with variable car- 
rier density has been derived. The phases present are the 
normal phase, where the order parameter is zero; the pseudo- 
gap phase where the absolute value of the order parameter 
is non-zero but its phase is random, and a superconduct- 
ing phase with a crossover quasi-2D Berezinskii-Kosterlitz- 
Thouless (BKT) region. The crossover region is bounded by 
the quasi-2D BKT temperature and the temperature for the 
onset of conventional long-range order (CLRO). The practi- 
cal observation of these regions however critically depends on 
the carrier density. At high densities both the pseudogap and 
BKT regions vanish asymptotically i.e. one obtains the stan- 
dard BCS picture. At intermediate densities the pseudogap 
phase is large but the BKT region negligible. Finally at very 
low densities both the pseudogap and BKT regions are size- 
able. An attempt is made to explain the behaviour observed 
in underdoped (intermediate densities) and optimally doped 
high-Tc superconducting compounds above their critical tem- 
perature. The transition to the pseudogap phase should also 
be regarded as a crossover. 
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I. INTRODUCTION 

It is widely accepted [fj]-[§ that the formation of the 
superconducting state in high-temperature superconduc- 
tors (HTSCs) is quite different from traditional su- 
perconductors. The strongly anisotropic, almost two- 
dimensional, structure of HTSCs compounds increases 
significantly the influence of fluctuations and the different 
types of order appear to be separated on different energy 
scales. One type of order may correspond to the forma- 
tion of a pseudogap (or lowered density of the states at 
the Fermi level) above the critical temperature T c of the 
superconducting transition. The anomalous behaviour of 
HTSCs ®M (including the behaviour of the spin suscep- 



tibility, resistivity, specific heat and photo-emission spec- 
tra) can then be interpreted in terms of the formation of 
this pseudogap j|||p|. 

There are many approaches for studying these ex- 
tremely complicated systems. In most models, whether 
they are discrete (Hubbard) or continuous, the systems 
are treated as two dimensional (2D). Even the simplest 
2D model with local nonretarded attraction between car- 
riers is still a very rich and complicated system. Its prop- 
erties beyond the mean field level are not well established, 
especially if one interested in a wide range of carrier den- 
sities and coupling strengths. 

Most theoretical approaches (see, for example, PJIc|l 
and the references therein) use the so called self- 
consistent T-matrix approximation, where all Green 
functions are "dressed". Although the self-consistency 
is an important feature [|To|— |T^] , we stress that even this 
approximation cannot adequately describe the formation 
of the Berezinskii-Kosterlitz-Thouless (BKT) phase and 
it is well known that this is the only type of order pos- 
sible in 2D (ll|. The features of the superconducting 
BKT transition were studied in jl6| although without 
self-consistency. 

It is important to note, however, that often real sys- 
tems (in particular HTSCs) are not strictly 2D systems 
and a pure 2D scenario cannot be applied to them di- 
rectly. The presence of a third spatial direction already 
permits the formation of a superconducting phase with 
long-range order (LRO). One cannot say apriori that the 
formation of this phase destroys the 2D physics abruptly. 
Nevertheless, using some estimations obtained more than 
20 years ago [17,18] for the quasi-lD metals, in principle 
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even for relatively modest anisotropies and carrier den- 
sities the standard BCS behaviour should be recovered. 
The low and intermediate density limit is not as clear 
and, given the relevance of this limit to HTSCs with their 
small and variable carrier density, it should be examined 
carefully. 

The purpose of this paper is to build the complete 
phase diagram of a quasi-2D metallic system with an 
arbitrary carrier density. We have considered both the 
conventional long-range order (CLRO) superconducting 
transition and the BKT transition. The accompanying 
publications were devoted to the calculation of the tem- 

tature for the BKT transition in the quasi-2D system 
, the temperature of the CLRO transition in the ex- 
treme Bose limit (see also review |p))). To explain 
the origin of the whole phase diagram and especially the 
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part related to the pseudogap region we recap here only 
the essential details from our previous works. The cal- 
culation for the temperature of the CLRO transition in 
the high density limit is new and we present here its full 
derivation (see also a preprint of authors devoted to the 
Coleman- Weinberg formalism in the theory of supercon- 
ductivity H). 

Indeed a quasi-2D BKT transition is only physically 
meaningful in the limit of weak three dimensionalization 
i.e. for a practically two dimensional system. In this limit 
we argue that for all carrier densities the BKT transition 
occurs before or at the CLRO phase. Certainly in the 
cuprates, there was until very recently little or no evi- 
dence for two superconducting phases and indeed it could 
be argued from experiment that there was only one su- 
perconducting transition. However in a recent paper ]22j | 
one indeed sees experimental evidence in underdoped Bi- 
cuprates for the onset of planar superconductivity before 
that of superconductivity perpendicular to the planes. 
This experiment which gives a temperature difference of 
the order of 10K thus supports the picture of the forma- 
tion of a quasi-2D BKT region prior to the formation of 
CLRO. 

We also want to argue that the calculation of the crit- 
ical temperatures for the quasi-2D BKT phase and the 
CLRO phase is very different both physically and math- 
ematically. The BKT crossover region has a 2D counter- 
part since in two dimensions Tbkt is finite; the CLRO 
phase does not since T c is zero p5| . Moreover at very low 
(physically unimportant) densities the transition temper- 
atures are straightforward to calculate and very different 
pof . The closeness of the two temperatures in the physi- 
cally important region (close to optimal doping) is a non- 
trivial result since the calculations are different in both 
the physical assumptions made and the mathematical de- 
tails of the approximations. We find it remarkable, given 
the highly anisotropic nature of the system, that the tran- 
sition to conventional LRO takes place at a temperature 
close to that of the BKT transition. 

We speculate that the quasi-2D BKT region corre- 
sponds to the crossover scale to three dimensional order. 
There may well be no sharp transition at the critical tem- 
perature for BKT phase formation. 

More importantly, the superconducting region is sep- 
arated from the normal phase by the pseudogap phase 
with nonzero neutral order parameter |l6| , |l9[ ]. Clearly 
the transition to this phase should be a crossover but our 
approximations are insufficiently accurate to give this be- 
haviour (see also discussion in jlj|). Furthermore this 
phase is two dimensional and large enough to explain 
the observed anomalies when the density of the carriers 
is relatively low (the underdoped case). Nonetheless the 
transition to the superconducting phase is almost directly 
into a state with the conventional three-dimensional long- 
range order at relevant doping levels. Thus we have a 
crossover with increasing temperature from three dimen- 
sional behaviour in the superconducting state to two di- 
mensional behaviour in the pseudogap phase similar to 



the picture recently proposed by p3| . 

The rest of the paper is organized as follows. In Sec. II 
we introduce the model and sketch the formalism ap- 
plied. This formalism is suitable for studying both the 
BKT and CLRO transitions. The equations for Tg KT , 
the temperature for the BKT transition in the quasi-2D 
system, have been derived previously [ jl9| and are simply 
sketched. The CLRO superconducting transition tem- 
perature for small and relatively large carrier densities 
is derived in Sec. III. In the high density limit we go be- 
yond the Nozieres and Schmitt-Rink approach |2^| which 
in the quasi-2D case, i.e. in the presence of a third di- 
rection, simply gives the BCS critical temperature. By 
employing the Coleman- Weinberg effective potential f23] 
we are able to calculate the Gaussian corrections to the 
critical temperature and we obtain the correct limiting 
behaviour T c — > as the system becomes increasingly 
two-dimensional. The full phase diagram of the system 
is discussed in Sec. IV and Sec. V presents our conclu- 
sions. 



II. MODEL AND FORMALISM 

Since the precise nature of the interplane tunneling 
in HTSCs is not yet known |2(J several different mod- 
els exist. Here we choose the simplest possible Hamilto- 
nian density which is widely used in the study of HTSCs 
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H= -V£(r) 



2mj_ m z d 2 
VV|(r)^(r)V|(r)V T (r), 



cos(zdV z ) + fi 



Mr) 
(2.1) 



where r = r±,r z (with being a 2D vector); Vv(r) is 
a fermion field, a =T,| is the spin variable; m± is the 
effective carrier mass in the planes (for example Cu02 
planes); m z is an effective mass in the z-dircction; d is 
the interlayer distance; V is an effective local attraction 
constant; /x is a chemical potential which fixes the carrier 
density nf \ and we take % = feg = 1 and m z » m±. 

In the weak coupling limit it is appropriate to replace 
the attraction constant V by the two-particle bound state 
energy in vacuum pq,E3| , 



£b = — 2Wexp 



4ird 
m±V 



(2.2) 



Here W is the width in the plane and the limit V — > 0, 
W — ► oo is to be understood. This replacement enables 
one to regularize the ultraviolet divergences in the theory. 
One can then define the dimensionless system parameter 



e = 



eg 

M 



(2.3) 



where e << 1 corresponds to Bose or local pair supercon- 
ductivity and e >> 1 to BCS superconductivity. Since 
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we have a quasi-2D system with a quadratic dispersion 
law in the planes the Fermi energy ep is given by 



e F = 



ttufcI 



(2.4) 



It was suggested in M that optimally doped HTSCs 
have e ~ 3 ■ 10 2 — 10 3 while conventional metallic super- 
conductors have at least e ~ 10 3 — 10 4 . 

The proposed Hamiltonian proves very convenient for 
the study of fluctuation stabilization by weak 3D one- 
particle inter-plane tunneling. The two particle (Joseph- 
son) tunneling has been omitted in (2.1) on the as- 



sumption that it is less important than the one-particle 
coherent tunneling already included. In certain situa- 
tions however Josephson tunneling is more important. 
In fact some authors consider the most important mecha- 
nism for HTSCs to be the incoherent inter-plane hopping 
(through, for instance, the impurity (localized) states or 
due to the assistance of phonons). Nonetheless we omit 
Josephson tunnelling here. However the single-particle 
dispersion relationship used here automatically incorpo- 
rates the layered structure of HTSCs which is a vital 
extension to the commonly used 2D models. 

It is significant that the large anisotropy in the conduc- 
tivity cannot be related to the anisotropy in the effective 
masses m z and m±. In particular, HTSCs with rather 
large anisotropy in the z-direction do not display conven- 
tional metallic behaviour at low temperatures [ pl| . How- 
ever this semiconducting behaviour is not directly related 
t o th e pseudogap phenomena ]3l| ] and the Hamiltonian 
( |2.l|) may be used to study the qualitative features of 
pseudogap opening. 

The standard Hubbard-Stratonovich method was used 
to study the Hamiltonian (pTl|). In this method the sta- 
tistical sum Z(v, fj, 7 T) (v is the volume of the system) is 
formally rewritten as a functional integral over the aux- 
iliary fields $ = Vtpiipi and <f>* = ^V'jV'j 

Z(v,n,T) = J V$V<f>*exp[-!3n(v,!i,T,$(x),$*(x))}, 

(2.5) 

where 

yj dr J dr\^{x)\ 2 -Tr LnG" 1 ^), $*(:r)] (2.6) 

is the effective action and x — r, r denotes both the imag- 
inary ti me t and the position r previously defined. The 
action (2J3) is expressed in terms of the Green function 
G which has the following operator form 



G _1 [$(x), $*(x)] = 



-Id r + T 3 



1 



2m± m z d 2 

-T+$(x) + T_$*(x). 



cos(idV z ) + fi 



where t 3 ,t± — (t% ± £t2 )/ 2 are Pauli matrices. Although 
the representation ( |2.5| ), ( |2.6|) is exact, in practical cal- 
culations it is necessary to restrict ourselves to some ap- 
proximation. The most convenient approximation for our 
purposes is the Coleman- Weinberg pj| (see also [ p2[ ) ef- 
fective potential in the one-loop approximation w hich we 
will use for studying LRO. The exact expression (2.5) is 
replaced by 

Z(v, fi, T, |$| 2 ) = exp[-(m pot (v, n, T, |$| 2 )], (2.8) 
where the effective thermodynamical potential 

n pot (v,[x, t, |$| 2 ) ~ n^ t ( v ,n, t, |$| 2 ) + nW(v,[x, t, |$| 2 ) 

(2.9) 

is expressed through the mean- field "tree-potential" 
n^ t (v, M , T, |$| 2 ) = Q(v, ft, T, *(x), **(a:))L *. = eoMt 



(2.10) 



and the one-loop (quantum) correction 



n^(v, fi,T, |$| 2 ) = -Tr Lnr; 1 + Ztt LnTI 1 (2.11) 



Here the Green functions T± are given by 

ri 1 (r,r) = 
(3S 2 Q 



± 



ps 2 n 



<5$*(r,r)5$(0,0) <5$*(-r,r)<5$*(0,0) 



<E> = <E>* = |<1?|= cons 



= v S(r)S(r) + tr[G(r, r)r+G(-r, -r)r_] 
±tr[G(T,r)r_G(-r,-r)r_]. 



(2.12) 



Strictly speaking one can only use this factorization into 
1 in the limit of small momenta and frequencies 



r+ x r 



(2.7) 



and hi gh c arrier density |33j] . While the representation 
fl2.5|), ( |2.6| ) is convenient for the studying of the LRO 
transition, the formation of the BKT and pseudogap 
phases demands a more subtle treatment. This is re- 
lated to the fact that these phases do not display LRO. 
This is evident in the original 2D model Jl6| where LRO 
is forbidden by the 2D theorems |H| . 

Thus one must avoid the presupposition about th e ex - 
istence of LRO used to write down the expression ([2.9]). 
The ill ega l step in the treatment based on representation 
(|2.5| ), ( pli| ) is the fixing of a definite value to the phase 9 of 
the complex orde r pa rameter $ which is performed when 
one obtains Eq. fl2.9|). To avoid this dangerous step one 
must use the modulus, p(x), and phase, 6(x), parameter- 
ization of the order parameter $(x) = p{x) cxp[i6(x)\ as 
first pointed out by Witten . This particular choice of 
parameterisation ensures that is single- valued with 
period 2ti. At the same time as this replacement by 
modulus-phase variables, one must reparameterize the 
Nambu spinor as ip a (x) — x<x(x) exp[i6(x)/2]. 
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As a result we have obtained 1 19 instead of ( |2.5| ), (|2 
the following representation 

Z{v,fx,T)= ( pVpD6exp{-(3n(v,n,T,p(x),d6(x))], 



(2.13) 

where the one-loop effective action ft(v, p, T, p(x), d6(x)) 
now depends on the modulus-phase variables and has 
been evaluated in Il9| where we have shown that 



n 



n kin (v, p, T, p, 09) + n^ t (v, n, T, p) (2.14) 



where the potential energy term in terms of p 2 is identi- 
cal to the mean-field potential in the BCS approximation 
but with |$| 2 replaced by p 2 . Thus T p , the temperature 
at which p = 0, is in this approximation identical to the 
BCS mean field temperature T C MF . However, although 
T p = T C MF in the mean- field approximation for p(x), the 
two temperatures have a very different basis, both math- 
ematically and physically p^Jl^ |. 

It has been shown that, in the model under consider- 
ation, the kinetic energy term reduces in lowest order to 
Ipjf the Hamiltonian for the classical spin quasi-2D XY 
model |35| (see also ]36|,[37]]) and we have simply used 
their expression for the BKT transition in the highly 
anisotropic case, when the vortex ring excitations are ir- 
relevant |35| . 

A self-consistent calculation of Tg KT as a function of 
rif also requires additional equations for p and p. The 
relevant equations are given by minimizing the potential 
VL p ot with respect to p and fixing the number density. 

The numerical investigation of Tg KT and T p has been 
carried out (see |l^|) and fortunately gives results not too 
different from the 2D case [[16i . 



III. SUPERCONDUCTING TRANSITION INTO 
THE PHASE WITH LONG-RANGE ORDER 

As discussed in Section II to study CLR O formation 
one should use the effective potential (2.9) obtained in 
the Gaussian approximation. 

There are two limiting cases where one can obtain an 
analytical solution, namely, the low and high density lim- 
its. For small concentrations of the carriers we have pre- 
viously obtained the critical temperature Q using the 
approach developed by Nozieres and Schmitt-Rink j24|] 
(see also |}8],[59|), while for the high density case we em- 
ploy a new method similar to that used in f37j| . 

The main point of is the solution of the system 
of n umb er and gap equations derived using the poten- 
tial (2.9) but with the one-loop corrections (2.11) taken 
on the critical line $ = <f>* = 0. This means that the 
effect of the fluctuations is only included through the 
number equation, while the gap equation is derived us- 
ing the mean field potential. Thus one always recovers 
the standard mean-field (BCS) gap equation 



1 

V 



rfk 



1 



(2tt) 3 2£(k) 



tanh 



2T r 



(3.1) 



where £(k) = k 2 _/2m — (m z d 2 ) 1 cosfc z d — p, while the 
corresponding number equation has the form 



n F (p, T c ) + 2n B (p, T c ) = n f . 



(3.2) 



One can see from (3.2) that the fermions are divided 
into two coexisting systems: fermi-particles, or unbound 
fermions with density np(p,T c ), and local pairs, or 
bosons with density ub{p, T c ). In the extreme Bose limit 
np = the temperature has already been shown to be 
the Bose condensation temperature of an ideal quasi-2D 
Bose-gas @|0) 

We now consider the high-density limit. This case is 
characterized by the condition that the Fermi surface is 
not disturbed by the attractive interaction, i.e. the con- 
tribution of the bosons to (3.2) is negligible and p ~ ep. 
Thus one need only study the gap equation. It is evi- 
dent that the Nozieres and Schmi tt-R ink approach which 
yields the mean- field equation (3.1) cannot describe a 
quasi-2D system adequately. Indeed it follows from the 
2D theorems that in the limit m z — > oo, when the sys- 
tem becomes two dimensional, the value of T c must go to 
zero. Clearly this can never be obtained from Eq. ( |3.1| ). 

To obtain the gap equation which does describe the 
quasi-2D system one must include the fluct uat ions. This 
can be done if uses the ef fective potential (2.S) with the 
one-loop correction (2.11) but without setting $ = $* = 



prior to taking the derivative with respect to (f>. 
Thus the gap equation takes the form 



^- F t ( W , M ,T c ,|$| 2 ) 



5|$| 2 

dnW( v ,p,T c ,m 2 ) 



0l#l 



(3.3) 



wh ere the quantum correction (2.11) which was omitted 
in ( |j.l| ) has now been included. 

There is a subtle p oint related to Eq.( |3.3| ). The value 
of T c defined by Eq.(3.3) should be less then the mean- 
field temperature, T^ b defined by Eq.([3.l|) (and given 
by flA.llf)) in the high-density limit p = e F > T C MF . 



However sketching the tree-potential OMf (T) as a func- 
tion of $ for T < T C MF yields the Mexican hat shape. 
This has a minimum at |<I>| 2 = $ 2 ni „ ^ (where $ m j n 
is simply the mean-field value for $) but a maximum at 
$ = $* = 0. Thus the role of the correction to the effec- 
tive potential in Eq.( |3.3[ ) is to transform the maximum 
of the tree-potential to a minimum of the full potential 

n pot m 2 ). 

Unfortunately, the one-loop correction f^ 1 ) defined by 
( p. 11 ) is ill-defined (complex) at the point of interest. 
The temperature at which the potential just becomes 
complex in fact gives the Thouless criteria of supercon- 
ductivity in BCS theory. However, in our treatment this 
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is only an indication that the one-loop approximation 
fails at the point <f> = <f>* = and is evidently related to 
the non-convexity of f^f (l^l 2 ) at this point. 

This situation is standard in quantum field theory if 
one considers the class of theories with tree-level symme- 
try breaking and an extensive literature exists |52|,fll| (l3) . 
In our case to describe a homogeneous state where |<&| 2 is 
uniform one should replace £l po t by the so called Gaussian 
PH or modified Q| effective potential which coincides 
with Opot in the region where the latter is well-defined. 

One could now find the modified potential, but we 
will use here, in our opinion, a more transparent con- 
sideration which will allow us to evaluate the one-loop 
correction not at $ = but in the region where it is 
well-defined and coincides with the modified potential. 
Our results can then be straightforwardly related to the 
modified potential of Weinberg E5] . 

Let us assume that T C <T C MF which means the point 
$ mi „ is close to zero. At this poi nt O ^f is surely con- 



vex and the one- loop correction (2.11) is real and well- 
defined. Thus for T C <T^ IF one can approximate Eq.( |3.3| ) 
by 



1 dnUf(v,e F ,T C: &) 



v <9$ 2 

1 dnW(v,e F ,T c ,<i> 2 ) 



<9$ 2 



(3.4) 



where we choose $ to be real and the value of $ m j n 
is simply the well-known mean-field BCS value f or <I> at 
temperature T (see e.g. 0] and Appendix A Eq. (|A.13|)) . 

Therefore to solve the approximated gap equation fl3.4| ) 
one has to calculate 



1 cmW(v,e F ,T,$ 2 ) 
v 



1 T 



JdKT ± (in n ,K] 



2(2tt)3^ ^ 



a$2 



(3.5) 



where K = (Kj_,if z ) and, starting from ( 2.12 ), one can 
obtain the Green's functions as a function of 4> in the 
momentum representation 

r i 1 (ffi„,K) = i + ^3 £ fdkx 

1— — OG 

[iwi - £_] [iuoi + ifi„ + £+] ± $ 2 



[w, a + e + * 2 ][(wi + n n ) a + ^ + * 2 ]' 



(3.6) 



where we have introduced the short-hand notation £± — 
£(k ± K/2) and Q„ = 2T:nT , uji = tt(2Z + 1)T are odd 
and even Matsubara frequencies, respectively. 

Since n 1 (0,O) = is simply the BCS gap equation 
it has solution <I> equal to the BCS value 3> OT j„ i.e. for 
$ = $mini ri 1 (0, K) has a zero at K = 0. This gives rise 



to the only singular term in ( p.5[ ) for $ = ^ m im namely 
the pole in r_(0, K) at K = 0. One can therefore write 



1 a0W(u,e F ,T,$ 2 ) 
v 



1 T 

2 (27r) 



^3 y dK r_(o,K)| $=$mir 



srz^o.K) 



a$2 



(3.7) 



In order to perform the calculations analytically we use 
the derivative expansion for the Green function TZ 1 , 



6[1 - cos^d]] 



(3.8) 



(see Appendix A, Eqs.(A.lO) and (A. 12)) and its deriva- 
tive 



arz 1 (Q,K) 



a$2 



mj_c 
2?rd 



(3.9) 



(see Appendix A, Eqs. (A.S) and (A. 4)) 



Substituting these two expressions into (3.7), one ar- 
rives at the following approximation 



1 9fiW(u,e F ,T,$ 2 ) 
1 T 



2 (27r) 3 



a$2 

dK- 



c 



iKi + 6[1 - cos K z d] 



(3.10) 



One can see that Eq.( 3.1(j| ) has no infrared divergencies 
due to the presence of the third direction (6^0). In two 
dimensions it would be infrared divergent as required by 
the 2D theorems jl5| . This equation also has an artificial 
ultraviolet divergence as a result of the replacement of 
the Green's function T" 1 by its derivative approximation. 
Thus one should introduce a rather natural ultraviolet 
cutoff (K™ ax ) 2 = 2m^$(T = 0) = 2m^2\e b \e F and 
integrate over the momentum K to obtain the expression 



1 dSlW(v,e F ,T,$ 2 ) 



where 



<9$2 



1 w 2 
AV2 el 



m± T 
2nd2e^ 



\1tlk\, (3.11) 



,d 2 ' 



(3.12) 



Substituting ( 3.11 ) into (3^4) one obtains the final tran- 
scendental equation for T c 



In 



T, 



T™ F 2e F 



ln/d = 0, 



(3.13) 



which may be rewritten in the following more convenient 
form 



T c = 2e F 



HT C /T t 



MF\ 



llnrel 



(3.14) 



Strictly speaking the equation ( |3.14 ) is only valid when 
T C <T^ IF . However, we intend to use Eq.(3.14) even 
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when T c < T^ IF , which will allow to qualitatively sketch 
the whole phase diag ram of the system. 

If one solves (3.14) numerically for reasonable width 
one obtains roughly linear behaviour T c ~ ep ~ nj for a 
wide range of intermediate densities (see below) . This de- 
pendence is observed experimentally Q for HTSCs sam- 
ples which have a Fermi surf ace. 

A second feature of ( |3.14| ) is that increasing \eb\ in- 
creases both T^ IF and T c . This is opposite to the be- 
haviour in the low density limit. Finally T c goes to zero 
as m z — > oo (w — > 0) as it must. This is very different to 
the limiting behaviour of Tg KT , namely Tg KT — > ?bkt 
in the limit m z — > oo. 

As stated abov e on e can also understand the approx- 
imation used in (3.4) in terms of the modified effective 
potential defined in Ji"^ . The modified effective poten- 
tial in 1 42 1 is defined as the minimum value for f2 given a 
homogeneous state where |$| 2 is uniform. The real part 
of this modified potential has the following form 

f2W(«, M ,T,|$| 2 ) = 

E / dK\nr ± \in n ,K), (3.15) 



2(2 



where the area T> of integration in the momentum space 
includes only positive mode s. One can see that ( 3.15 ) 
indeed coincides with (2.11) when f2 W is well-defined. 
Furthermore t he m odified potential ( 3.15 ) leads to the 
gap equation (3.4) which was considered above as the 
approximated one. 

In the region $ < Q m in the modified effective potential 
considered above differs from the traditional effective po- 
tential, Q e ff(Jf>), which is defined as the minimum value 
for SI such that the space average of Q(x) is given by $. 
It can be shown that the conventional effective potential 
is in fact the convex envelope of the modified effective 
potential and is real and convex everywhere. However 
for $ < & m in it describes an inhomogeneous mixed state 
where the value of $(x) is not uniform in space. One can 
readily understand that the modified and not the original 
potential is relevant for the superconducting state. 

There is, however, the difference between our interpre- 
tation of the modified potential and t hat of fl42[ . In 
the homogeneous state described by ( fr.lE ) is considered 
as decaying and the rate of the decay is related t o ncg - 
ative modes of ( 2.11 ) which are not included in ( 3.15 ). 
It is physically obvious that there is no real decay of the 
homogeneous superconducting state with $ < $ m i„ for 
T < T c although we have not been able to prove this rig- 
orously. The absence of decay is in agreement with the 
interpretation of Q] although it should be stressed that 
the modified potential discussed here is not identical to 
the Gaussian effective potential in El| . 



IV. THE PHASE DIAGRAM OF THE SYSTEM 

We have now obtained solutions for T Cl Tg KT and T p 
and can build the phase diagram for the quasi-2D system. 
This diagram comprises the normal phase, the pseudogap 
phase, the superconducting BKT region, and the CLRO 
phase with superconductivity in the bulk in order of de- 
creasing temperature. The practical observation of these 
regions (i.e. the temperature range for which they are 
present) is however critically dependent on the carrier 
density. 

i) We first reiterate results at very low carrier densities 
e <C 1 i.e. for a Bose liquid of isolated local pairs. 

Firstly, the critical temperature is linear in ep T c *~ ep 
(or T c ~ rij as expected in the 2D case) p0[ . For large 
anisotropy the proportionality constant depends only on 
the particle density and the width in the z-direction. This 
constant is less than 1/8 for half-bandwidth 



w < 



2\/2exp2 2y/2 exp i 



(4.1) 



where $(T = 0) is the zero-temperature energy gap. Re- 
call that in the 3D case T c ~ nj 3 |38| . Secondly, contrary 
to the case of a 3D superconductor where T c is indepen- 
dent of V @ , in a quasi-2D system T c does depend on 
V. In particular T c decreases with the growth of V . The 
reason for this is the dynamical increasing of the compos- 
ite boson mass in the third direction. Thus, the growth of 
\eb\ (or equivalcntly of V) "makes" the system more and 
more two-dimensional in this model of a quasi-2D metal 
with a local four-fermion interaction. It is interesting to 
note that a decreasing T c can also take place in the case 
when the local pairs are bipolarons p5| . 

Furthermore the temperature Tg KT always lies above 
T c in this limit for reasonable bandwidths (see Eq. (4.1) 
and strongly anisotropic systems Pq,f9| since 



^BKT > ^BKT 



(4.2) 



In addition the BKT region and LRO order phase are 
of roughly equal size so one expects to see an extensive 
region where the interlayer tunneling is insufficient to 
produce CLRO. 

This range of carrier densities is not however experi- 
mentally important because the anomalous behaviour of 
HTSCs is observed when the Fermi surface is still present. 

it) We next present results for high carrier densities 
(e~> 10 3 ). 

At high carrier densities T c approaches T C MF asymp- 
totically so that the condition T C <T C MF is satisfied and 
equation (3.14) for the critical temperature is indeed 
valid. F urthe rmore one can expand the logarithm in 
equation (3.14) to give 



T n =T n \l — 



T p \ \iik\ 
2e F 



2>KT 



(4.3) 



G 



where T C MF has been replaced by the temperature T p 
(which is identical in the approximation used in this pa- 
per). 

For these carrier densities one also finds that 
p(2bkt)/Tbkt <C 1 which yields the following estimation 
for the BKT transition temperature 



T p > Tg KT > Tbkt — T p ( 1 — 



4Tp 



e>10 3 . (4.4) 



From (gj) and (|J) it follows that Tg KT > T r for rea- 
sonable anisotropies, bandwidths satisfying (4.1) and for 
particle density, e > l/2(e/2) 8 which is clearly satisfied in 
the high-density region. Furthermore both T c and Tg KT 
tend asymptotically to T p — T C MF in the high density 
limit. Thus both the pseudogap and BKT regions vanish 
asymptotically in this limit and one recovers the standard 
BCS behaviour with a transition directly into a supercon- 
ducting state with CLRO at the BCS critical tempera- 
ture as obtained by |17| in the quasi- ID case and in the 
weak-coupling limit in 2D ||[l6 46 1. 

m) The physical interest lies at intermediate densi- 
ties (e ~ 10 — 10 3 ). Optimal doping corresponds to 
e ^ 3 • 10 2 — 10 3 |32] an d we show this range of den- 
sities in Fig. 1. The precise value of e for optimal doping 
is however strongly dependent on the anisotropy of the 
system. For these relatively high densities the condi- 
tion T^ KT , T C <T C MF is still satisfied and equation ( ]3.14| ) 
is a good approximation. For example, for e = 700, 
(T p — T c )/T c sa 0.15, indicating that both the pseudo- 
gap and BKT regions are small. One can see from Eq. 
(O) that 



Tp — T c 



0. 



(4.5) 



This gives an asymptotic disappearance of the pseudo- 
gap since the fraction of the temperature range in the 
pseudogap phase goes to zero. This corresponds to the 
experimental observations. Furthermore one finds that 
T c is very close to Tg KT which is in agreement with 
the experimental observation of at most a very narrow 
crossover region in temperature from two-dimensional to 
three-dimensional superconductivity reported in [ ^7|j22| 
for optimally doped cuprates. 

Due to the fact that the Fermi surface is still present 
for the densities where the pseudogap develops, one 
might expect the pseudogap pheno mena to occur at 
e ~ 10 — 10 2 . Extrapolating ( 3.14 ) to these densities 
one finds that T c is always less then Tg KT . More essen- 
tially {Tp — T^ KT )/T p ~ 1, i.e. the pseudogap phase is big 
enough to make it observed and sufficient to explain the 
observed anomalies in the underdoped compounds. On 
the other hand (T^ KT - T c )/T c < 1 so that the BKT re- 
gion again remains relatively small. In this density range 
these two temperatures are so close that one cannot see 
the difference on the phase diagram. The phase diagram 
is thus unchanged from that shown in [ |l9| ] since the lines 
for T c and Tg KT coincide. Thus for all carrier densities 



appropriate to HTSCs the superconducting transition is 
practically directly to the phase with CLRO. Since the 
temperature range for the BKT region is only a few de- 
grees it has until recently [E2| lain within the observed 
critical region. 

We note that our results for T c imply that the quasi-2D 
BKT region is decreasing in size with decreasing density. 
We believe that this is not the case physically particu- 
larly since it has been shown that in the Bose limit the 
quasi-2D BKT region is large. Thus we believe that this 
is an artifact of our approximation since it becomes less 
accurate as T c moves away from T p which occurs with de- 



creasing density (see discussion after Eq. (3.14)). Thus, 
although one can say that this transitional temperature 
range is small, our approximation is not adequate to de- 
scribe its dependence on the carrier density. 



V. CONCLUDING REMARKS 

We have derived a phase diagram for a superconduct- 
ing quasi-2D system in which one has four regions. These 
are the normal phase, the pseudogap region, the BKT 
region and the phase with conventional bulk supercon- 
ductivity (CLRO). 

For all realistic carrier densities we find the BKT re- 
gion to be only a few degrees. This implies that it lies 
within the critical region and is difficult to observe experi- 
mentally. In fact Tg KT is slightly overestimated (since for 
example the effect of vortex rings is omitted) . In contrast 
the value of T c is underestimated by the approximation 
used since only one loop corrections have been included 
p7[ . Thus the BKT region, if present, is even smaller 
than calculated. 

In contrast the pseudogap region is substantial for 
the underdoped region where the unusual superconduct- 
ing condensate formation is responsible for the observed 
anomalies in underdoped HTSCs. At optimal doping the 
pseudogap phase is also negligible and one recovers the 
standard BCS behaviour. 

With the calculation of T c we are able to present a pic- 
ture in which two dimensional behaviour predominates 
in the pseudogap phase but where the superconducting 
phase is three dimensional in line with the recent result 
of & 
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APPENDIX A: THE CALCULATION OF THE 
GREEN FUNCTIONS T± 

We derive here the Green functions r±(if2„, K) in the 
derivative approximation for the la rge d ensity (/i ~ 3> 
T) limit starting from expression Since $/T C 1 

for T<T C one can use the "high-temperature" approxi- 
mation 



1 



1 



$2 



C 2 + $ 2 



which gives 



1 



r ± 1 ( i r»„,K)-- + r £ 

I— — oo 



1 



(2^)3 X 

2 $2 ± $ 2 



(A.l) 



(A.2) 



where to lowest order one can set f2„ = K = in the $ 2 
term. Performing the summation over fermion Matsub- 
ara frequencies one arrives at the following result 

r^ 1 (in n ,K) = r- 1 ^, k) + ^c(2$ 2 ± $ 2 ), (a.3) 



where in for large densities 



8tt 2 T 2 



(A.4) 



and T(iCl n , K) is the Green function of the order param- 
eter fluctuations in the normal ($, $* = 0) state: 



V ' V 2 J (2irf 



tanh -i + tanh — - 
2T 2T 



(A.5) 



The 2D limit (to = 0) of the Green function was studied 
by Shovkovy et al. jl8],fl9| and we need to generalize the 
expression to the quasi-2D case. Substituting the cou- 
pling const ant V expressed via the bound state energy 
(2.2) into (A.5) and performing simple transformations 
on the hyperbolic functions one arrives at 

And \e b \ 



1 



r 2W r-2-K 

i dx / dt 

Sn^d ./n .In X + S(t, K_l, K z ) - lU n 



x tanh 
+ %nH. 



o J o 

x + s(t,K ± ,K z ) 
AT 

I-2W n2-n 2 

dx / dt— 







(I 



7T 



tt/2 



dtp 



x + s(t, K±,K Z ) - 



, ,x + s(t,K±,K z ) 

x tanh —— x 

AT 



cosh ^^^' COS ^ -l 



cosh x + a(t,K ± ,K z ) + cQsh d(t,K ± ,K f ,co8<p) 
where we introduced the short-hand notations 



, (A.6) 



s(t,K ± ,K z ) 



K 2 , „ (K r d 
2w cos t cos 



4mj_ 



2/i 



and 



„ , k±K± cos as „ . . ( K 7 d 
d[t, K± , A z , cos </j) = h 2w sin i sin 



Using the identities: 
1 



1 



x + s(t,K x ,K z )-i^ n x + s(t,K x ,K z ) 

+ [x + s(t, K ± ,K z )][x + s(t, K ± ,K Z ) - ifi„] 
and 
ra dx 



(A.7) 



4a7 



tanh a; = In 

X 7T 



dx 



(tanh x — 1) , at a > . 

(A.8) 



(In 7 i s the Euler constant), one can take the limit W 
oo in (A.6) and come to the expression: 



And |e 6 |7 



8ir 2 d Jo 



dt 



dx 



• tanh 



o 



X 



s(t,K ± ,K z 
AT 



m_i_ 



2tt 



eft / dx 
<< Jo 

1 



[a; + s(t, K±,K z )][x + s(t, K±,K Z ) - iO n ] 



x tanh 



AT 

2-rr 

dt 

ii Jo 



dx — 

7T 



x + s(i, A'x, if 2 ) - ift n 



x + s{t,K ± ,K z ) 

x tanh x 

AT 



cosh d(tJC H K Zl co^_ 1 



cosh x + s{t,K u K z ) + cQsh d(t,K ± ,K r coa<p) 



(A.9) 



One may e asily show that in the 2D limit w — > the 
expressio n (IA.9 T) transforms into the Green fu ncti on ob- 
tained in ||48||. For example, the last term of ( A.9 ) gives 
the exact representation for the imaginary part of this 
function. 
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We, however, are interested in the derivative expansion 
of r (0, K) for the anisotropies w/T <C 1 and large 
densities which has the following form 



r- 1 (o,K) = 



2?rd 



In 



T 



aKi 



b\l 



where 



T. 



MF 



: y/2\e b \e F 



cos K z d] 
(A.10) 

(A.ll) 



is the mean-field transition t em perature obtained from 
the mean- field gap equation (fo.lj); 



7C(3) e F 



b = 



7C(3) w 2 



(A.12) 



(47r) 2 m ± T 2 ' (An) 2 T 2 ' 

Again the expression ( A.10| ) is in the correspondence with 
the Green function used in [^9|. We note also that the 
anisotropy a(e F ) = bd 2 /2a is equal to the anisotropy that 
was found in the quasi-2D BKT model . 

Fi nally we sh ould s tress that if one substitutes into 
(A.3) (see also ( A.10 ) the well-known mean- field BCS 
dependence (see e.g. ]44|) 



8tt 2 {T c mf ) 2 
7C(3) 



$ 2 . (T) 

III ! II \ J 



1 



t mf) 



(A.13) 



of th e gap on the temperature near ', one obtains 
(U) which shows that T + (0, K) > 0, while T_ (0, K) < 
with the pole at K = 0. 
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FIG. 1. T c , T,| KT and T p = T C A/F versus noninteracting 
fermion density. The regions of the normal phase (NP), pseu- 
dogap phase (PP), BKT phase and long-range order (LRO) 
phase are indicated. We assumed that m z /mj_ — 100 and 
(m/hl)-^!. 
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